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■ Abstract 

We investigated the non-leptonic weak interaction in magnetic field. We dis- 

^ I cussed an improvement of previous method to analytical work out the rate for weak 

Q\ [ field case.Our result easily goes over to field-free limit.Then we calculated the re- 

. action rate in strong magnetic field where the charged particles are confined to the 

^ ' lowest Landau level.A strong magnetic field strongly suppressed the rate,which will 

O ■ be foreseen to affect viscous dynamics in SQM .We also derived a few approximation 

o 



formulae under given conditions that can be conveniently applied. 
PACS numbers: 97.60.Jd, 12.38.Mh, 95.30.Cq, 98.35.Eg 

1 Introduction 



^ I The composition of comparable number of u,d,and s-quarks are known as strange quark 

matter(SQM) would be stable or metastable configuration of hadronic matter. The bulk 
SQM is the /^-equilibrium system determined by such a serious of weak process: u + d ^ 
s + u, d ^ u + e + Ve ,s ^ u + e~ + Ve [IJ. There has been a lot of interest in the 
study of the reaction rate and their astrophysical relevance|l2l|3l|4l|5l|6l|7|. In the interior 
of neutron stars,the semi-leptonic reaction is devoted to the cooling of stars while the 
non-leptonic one to the damping of instability of rapidly rotating stars. We here concern 
about the non-leptonic process. The previous work|I^|3 were carried out under the usual 
situation of zero magnetic field. However, the strength of the surface magnetic field of 
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a pulsar is typically of order lO^^G which concluded from the observation data. Some 
magnetars are observed to have magnetic field of 10^^ ~ lO^^G . Considering the flux 
throughout the stars as a simple trapped primordial flux, the internal magnetic field may 
go up to lO^^G or even more IS 13. In spite of the fact that we do not know yet any 
appropriate mechanism to produce more intense field, the scalar viral theorem indeed 
allows the field magnitude to be as large as lO^^G. Therefore,it is advisable to study the 
effect of the magnetic field on strange quark matter including the calculation of the rate of 
the reactions. The quark Urea process have been discussed in a few of works (HI |3- Here 
we focus on the non-leptonic process leading to bulk viscosity in strange quark matter. 

u{l) + d ^ u{2) + s (1) 

In magnetic fields, a classical charged particle will have a transverse cycle motion. In 
quantum mechanism, the transverse motion is quantized into Landau levels. The quanti- 
zation effect is important when the magnetic strength is equal or larger than the critical 
value -B^-' = ni^c' /{qih) defined by equating the cyclotron energy qB/{mc) to mc^ . 
where rrii and qi denoted the mass and charge(absolute value)of the particle ifTOlfTll . fi,k 
and c denotes the Planck constant and Boltzman constant and velocity of light respec- 
tively, which are taken to be units below. We have considered a wide range of magnetic 
fields in our study, from "low" to "very high" magnetic fields(i? > 10^^(j)called respec- 
tively "weak" and "strong" field. Since u-quark mass in magnitude of order is smaller 
than d-quark and s-quark,we assume quantization effect on u-quark is important but that 
of other flavors are negligible[4|. In the weak field strength situation,we deal with the 
calculation of the rate through some approximations: (i)the unaffected matrix element; 
(ii)the approximating free-particle motion direction. In degeneracy,the strong magnetic 
fields,i.e. 2qB > pp^, force u-quark to occupy the lowest Landau ground state, where pp. 
is the Fermi momentum,we should use the exact solution of Dirac equation in magnetic 
field to evaluate the reaction rate. 

This paper is organized as follows. We solve the the dirac equation in part2. We 
consider the case of weak field in section 3. And then give the reaction rate in strong field 
in section 4. Finally, the results are discussed and a summery is made in section 5. 
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2 The solution of Dirac function in magnetic field 



We first solve Dirac equation in the presence of a magnetic field B.Let the uniform mag- 
netic field B along the z-axis,and we choose the asymmetric Landau gauge 



A = (0, Bx, 0) 



(2) 



so that the four-dimensional polarized wave function can be expressed in term of station- 
ary states in the normalization volume V — Lj-LyL^. In the magnetic field, the u-quark 
wave function for fermions ultra-relativistically reads 



exp[—iEt + ipyH + ipzz] 
^2E{E + m)LyL, 







\^ -i^2uqBI^_i.p^{x) J 



(3) 



for spin up, and 



■(/)_(t,x) 



exp[—iEt + ipy + ipz] 
^2E{E + m)LyL, 



\ 



{E + m)I,^i,p^{x) 
iy/2iyqBI„.p^{x) 



(4) 



for spin down cases, where the energy E — ^Jpl + + 2vqB , and v is the Landau 
level. 

It is degenerate and can be expressed as other quantum numbers such as, I is the orbital 
quantum number and s is the spin quantum number. 



(5) 



q is the charge of u-quark. In Eqn{3), (4), 



and 



l^qB 

where H^, is Hermite polynomial 



(6) 



(7) 
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3 The rate in weak magnetic field 



If the matrix element for the reaction(l) is determined, we can express the rate per volume 
ofreaction(l) asElll: 

T{md su,) = y [n / jSf^}\^s?S{27:Y6\P, + P, - ~ P.) (8) 

where the phase space integrals are to be calculated over all particle states, the statistical 
distribution function S = — /2) (1 — /«), the quark in equation(l) are described by 

Fermi-Dirac distributions in the form 

f,{E,) = [l + exp{^^j^)]'\z = 1,2, d,s (9) 

where /ij are the chemical potential. 

We consider the reaction (1) when the magnetic field is not strong enough to force the 
quark into the lowest Landau level. Previous study iQ shows that the matrix element for 
the weak process remains unaffected and only the phase factor is modified. 

The matrix element summed over final spins and averaged over initial spins is given 
bylH: 

|M,|2 = QAGlsm^e,cos^e,{Pi ■ Pa){P2 ■ P,) (10) 

here Pj = {Ei — pj is the four-momentum of the quark i and Gp = 1.166 x lO^^^Mef ~^ 
is the Fermi constant , 9c is the Cabibbo angle(cos^ 9c = 0.948).We neglect the masses 
of up and down quarks, then Ei = pi,E2 = P2, Ed = pd, Eg = {pi + mlY^'^. so the 
four-momentum products can be written as||2l|3l 

(Pi ■ Prf)(P2 ■ P,) = E^E^EdEsil - cosew)(l - If cos^^2s) (11) 

where 9ij denotes the angle between quarks i and j . 

Now we can calculate the rate of weak process in the weak magnetic field by replacing 
the u-quark phase space factorQ 

so the reaction rate can be write as: 

T{uid ^ SU2) = -—Glsm'9cCos^9c{eBy ^ (2 - S,,J J] (2 - 5,, J 

J Pddpdpldps j dpu J dp2zS6{Ei + E^ - E2 - Es)I (13) 



where 

d,s 

1=1 {l[)dQ,{l - cosei,)(l - ^ cos e2s)6\p, + p, - P2 - pJ (14) 

i ^ 

As we known ,the angle between two vector is a function of the respective inclinations 
and azimuth angle in spherical coordinates, ie. 

cos 9id = cos 9i cos 9d + sin 9i sin 9d cos(v9i — ipd) 

cos 6^28 = cos 6^2 cos 9s + sin 6^2 sin 9s cos(v92 — fs) (15) 

In magnetic field,the angles of the polarized quarks,6'i, 92, (fi, (p2 vary with Landau level 
(z/i, z/2) and z-component momentum {piz,P2z) ■ Therefore the integrals and summations 
in Eqn{12) become very difficult due to the coupling of variables u and pz into the inte- 
grated function. We need to make a improvement on chakrabaty's approach. Fortunately, 
we can approximately regard the angles as independent variables for the weak field sit- 
uation and then the integrals and summations decouples each other, because the kinetic 
quark direction should have only a small deviation from the field-free case,although the 
modification of the absolute value of the momentum considered under the situations.lt 
completely coincides with free-particle's matrix approximation described by Eqn{lQ). 
We immediately have: 

where Y,' and /' denote the summation and integral independent of the direction angles. In 
this approximation , which may be called free-particle direction average(FDA). Eqn(13) 
and Eqn(14) read: 

T{u^d ^ SU2) = jttteGI sin' 9, cos' 9,{eBf J] (2 - <5,,,J ^ (2 - J 

l^^/ ui=0 U2=0 

J Pldpdpldps J dpi, J dp2zS6{Ei + Ed-E2- Es)I' (17) 

where 

J' = y ( n )c?n,(l-COsM(l-^COS^2s)5'(Pl+Prf-P2-pJ (18) 

Since fXi >> T ,only those fermions whose momenta lie close to their respect Fermi 
surface can take part in the reaction. We use the method in ||3|7|to complete partially the 
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integral of /' through 



/u X 
^^exp(zp-x) (19) 



can be written as: 

2^7r2 f dx 

/' = —['[[sm{pF^x) + asm{pFiX)sm{pF^x)f{pF^x)f{pF^x) + 



sm{pF^x)sin{j)F,x)f{j)F^x)f{pF^x) + a /(m^))] (20) 



Pf^Pf2PFs 



-J (21) 

where a = f{pFiX) = cos{pFiX) — "^^^J"'^^ and the integral J is defined through 
Eqn{20) and ii^g?2(21),can be calculated numerically. 
The net rate of transforming d-quark into s-quark is |I3 

r{d ^s) = [l- exp(^^^^)]r(nirf ^ SU2) (22) 

Using the method in |I71[T2| and substitute Eqn{17), (18) and (21) into (22) we can get 

T{d ^s) = sin2 9, cos2 e^igBf {2 - S,^,o) E"'(2 - S,2,o) 

ul=0 u=0 

^^^^ =A^(A/x2 + 47r2T2)J (23) 



//if - 2uiqB^fii - 2u2qB 
where 

iyimax = Int{—^),i = 1,2, Aij = i^d- l^s (24) 
2gi3 

when the -B ^ the sum can be replaced by integral of v and then the result goes over to 
the expression: 

r(d^s) = ^Gl sin^ cos^ OcfilfilfisAfiiAfi^ + 4n^T^)J (25) 
It is just the field-free case. 



4 The rate of weak process in strong magnetic field 

Now we consider strong magnetic field effect on the non-leptonic weak interaction in this 
section.In the case of strong magnetic field that B > B^\ all u-quarks occupy the lowest 
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Landau ground state with the u-quark spins pointing in the direction of the magnetic 
field. We treat other flavors in this process as free particle which do not affected by the 
magnetic field. The matrix element for the reaction reads: 

G F sin 9 c cos 9 c 



^/2 



d,s quark treat as free particle ,then we get 



(26) 



A = —exp(-iPi ■ r)Ui 



( 1 \ 



lEi + rrii 



2E, 



Ei+rrii 
Ei+nii 



J 



V'2 

where 



LyLz 



TT 



Uo 



pE2{E2 + mu,) 



( Ei + m\ 


V J 

( E2 + m\ 


V J 



(27) 



(28) 



where Pi denotes four-dimensional momentum,^ = d,s. 

Consider v = 0,the wave function Eqn{3){A) of u-quark can be writen as 

^1 = — exp{-iEit + ipy^yi + ipz^zi){—)^exp[-^{-x + ^fpi (29) 
^LyL, TT 2 qB 



= J— exp{-iE2t + ipy,y2 + ip,,Z2){^y^exp[-^{-x + ^f]U2 (30) 



(31) 



(32) 



for spin up. 

We now calculate the matrix elements squared and summed over the initial state and 

average over the final state.We use the Eqn{29), Eqn{32) and Eqn{33) to express corre- 
sponding matrix element: 
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|M,|2 = [f727^(l - 75)[/,[7,7'^(l - 75)t/i][f727M(l - 75)t/dt/s7'^(l - 75)^^1]^ 
And then imediately get: 

|M,|2 = -^—-^{E,-p,,)\E2 + p.,)\E, + pJ{E,-pJ (33) 

Here we set = m„ = ,and then p^^ < 0,Pz.^ > O,otherwise jM^p = 0. Carrying out 
the integral ,we obtain 

5 {E2 +Es-El- Ed)S{Py^ + Py^ - Py^ " Py,)5{p,^ + " " ) (34) 

The rate per volume of the reaction is given by: 

r(„,<i ^ = / i:% hiij,^^ ^irfp„ 

/CO r /-oo r 

7^dp-z. / ^cip.,|M|Vi/d(l - /2)(1 - fs) (35) 

where the factor rii = Ud = Q comes from 2 spins and 3 colors. Since only left-hand 
helicity states of ui — quark couples to the W~ (W^ is the mediate of the reaction), a 
factor of I is shown in the Eqn(35). The integrals over dpy-^ and dpy^ can be carried out 
using the y-component momentum delta function. The integrals over dp^-^ and dp^^ are 
converted into dEi and dE2 respectively HQ. 

Since fXg >> T,only those momenta lie close to their respective Fermi surfaces can 
take part in the reaction. As the z-component momentum conservation,we can get 

PFi + PFi COS 9d - Pf2 - Pf, cos 6s = (36) 

We approximately set pp^ = Pf^ = PFs near the equilibrium.Then we can get cos6'd — 
cos Og = 2, so we can take out the integral over d-quark and s-quark's momentum space. 
Substituting Eqn(35) into Eqn(22), we obtain 

, ^ _ Gj, sin^ ^ecos^ e,{qB) ^ (2/i^)^ - {jJd + PFj^ 

(3/i, + 6/i„ - fid)filAfi{Afi' + in'T^) (37) 



where A/i = //^ - ij,s,Pf, = \I ^xl - ml. 



5 Discussion and Conclusion 



We have done calculations of the rate of non-leptonic quark weak process in magnetic 
field and give the analytic solutions under weak-field and strong-field approximation. 
Based on these results , we express the net rate of the transforming d-quark to s-quark 
in a unified form: 

r{d ^s) = Tkirib, g5)A^(A/x2 + 4tt^T^) (38) 

where k takes 0, L, H denotes zero-field, weak-field and strong-field cases respectively.In 
according with the formula(25), the result go over to the field-free case when the magnetic 
field strength vanishes. We thus have: 

ToK) = ^Glsin^^ecos^^el-)^ (39) 
meanwhile for 2quB < Ai^,the Tl reads simply from Eqn(23) 



144 ^9.9^ , ^ „ , vA 



TLirib, qB) = —G^ sin^ ^cCos^ OcQuBnbJ — j ■ r (40) 



rn 



where J as a function of fi define in figure(l).We find Tl has a small deviation from 
field-free case.The comparison is made in figure(2). 

The simplified formula of Th {rib, qB) need slightly complicated calculations for 2quB > 

In the referred SQM (//j ^ mi,only u-quark polarized)here , the number density of 
the components read: 

nd,s = ^ (41) 

^ ^^^^ 

/^d = /^s = = - /^e (43) 



the /3-equilibrium 
the charge neutrality 



2n„ - Ud - ris - Sue = (44) 
the baryon number density conservation 

"-6 = ^ {rid + nu + ris) (45) 



should be satisfied. Then we solve these equation numerically to obtain the chemical 
potential of the quarks with respective magnetic field. 

Figure(3) shows that the chemical potentials of quarks are nearly equal under this 
situation. Combines £'5n(41)(42)(44) and (45) in the approximation of ~ = A*s = 
fjL ,11 can be solved analytically through the algebraic equation 

//^ + -quBix - —Ub = (46) 

and then r//(nb, qB) is arrived at 

Gl sin^ Oc cos^ e^^qB) q^B 



(-en^TT^ + ^36n^7r4 + (g„S)3)l]3 (47) 

Figure(4) give a comparison of Eqn{A!7) with Eqn{37). They fit well each other with 
a small error. 

The formula (47) will reduce to 

= quBub (48) 

for 2qB ~ /x^,and tend to the limit 

IGTrGl, sin^ 0^ cos^ 
^ ^ ~ ? — m2 

{quBY 

for 2qB > //^ 

Figure(5) give a comparison of Eqn{Al) with Eqn{48) and £"571(49). We find they 
will be some good approximations in future realistic application. 

Under the consideration that u-quark is polarized but effect of other flavors are neg- 
ligible, we investigate the influence of magnetic field on the non-leptonic rate although 
the result for the weak field case has a small deviation from the field-free case ,we give 
an analytical treatment of the weak reaction which can be extended to the calculations 
of other reaction process. However, the strong magnetic field can extremely suppress the 
rate.It is possible to lead the decrease of bulk viscosity in magnetized SQM.This may 
have serious implications for compact star and pulsar dynamics. In fact, we should also 
distinguish d-quark and s-quark in calculations to obtain refining result for applications. 
It is our future works. 

This work is supported by NFSC under Grant Nos.90303007 and 10373007,and the 
Ministry of Education of China with project No.704035. 
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Figure 1 : numerical result of J , given by Eqn{20) and (21), show as as a function of for 
various value of the parameters m^. The solid curve is for = 80Mev,The dot curve is 
for nis = lOOMe^;. 
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Figure 2: Derivation of from Fq with B when rii, — 0.2 fm 
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Figure 3: numerical result of // as a function of B which get from £'gfn(46)(47)(48).The 
dot curve is for /x„ when = 0.2/m~^.The solid curve is for the case when /j, — — 
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Figure 4: Th as a function of B when rib — 0.2 fm ^.The solid curve is for the result of 
Eqn{37) and the dot curve is for the result of Eqn{Al). 
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Figure 5: F/^ as a function of B when Uh = 0.2/m^^.The solid curve is for the result of 
Eqn{A7) and the dot curve is for the result of Eqn{A8),the dash dot curve is for the result 
of Eqn{49). 
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